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(N ■ ABSTRACT 

Radiative torques, due to the absorption and scattering of starlight, are thought to 
play a major role in the alignment of grains with the interstellar magnetic field. The 
absorption of radiation also gives rise to recoil torques, associated with the photoelectric 
effect and photodesorption. The recoil torques are much more difficult to model and 
Q_|. compute than the direct radiative torque. Here, we consider the relatively simple case 

q ! of a spheroidal grain. Given our best estimates for the photoelectric yield and other 

£j \ relevant grain physical properties, we find that the recoil torques contribute at the 

~ 10% level or less compared with the direct radiative torque. We recommend that 
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the recoil torques not be included in models of radiation-driven grain alignment at this 
time. However, additional experimental characterization of the surface properties and 
photoelectric yield for sub-micron grains is needed to better quantify the magnitude of 



1. Introduction 



. these torques. 

, Subject headings: ISM: dust 
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As starlight traverses the dusty interstellar medium (ISM), it acquires a partial linear polar- 
ization, due to dichroic extinction. Grains are non-spherical and are more effective at attenuating 
radiation polarized along their "long" axes, which tend to align along the same direction in space. 
Despite decades of investigation, the alignment mechanism is still not fully understood; see Lazarian 
(2003), Roberge (2004), and Lazarian (2007) for recent reviews. 

Radiative torques, due to the absorption and scattering of starlight, appear to play a major role 
in grain alignment (Draine &: Weingartner 1996, 1997). If the incident radiation field is anisotropic, 
then radiative torques can directly align grains, at least for the few grain shapes that have been 
investigated so far. A detailed investigation, incorporating a large suite of shapes, is needed in 
order to fully assess the viabilty of radiation-driven alignment. Although the alignment is not due 
to magnetic torques in this scenario, the observed polarization is nevertheless correlated with the 
interstellar magnetic field direction, since the grains have magnetic dipole moments lying along their 
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spin axes and thus precess rapidly about the field direction (Martin 1971; Dolginov & Mytrophanov 
1976). 

In addition to the direct torque due to absorption and scattering, there are also recoil torques, 
associated with photoelectric emission and photodesorption. In the latter process, adsorbed atoms 
or molecules (i.e., gas-phase particles that have stuck to the grain surface after colliding with it) are 
ejected back into the gas after an absorbed photon breaks the bond between the adsorbed particle 
and the surface. So far, photoelectric and photodesorption torques associated with anisotropic radi- 
ation have not been considered. However, one expects that they could be significant compared with 
the direct radiative torque, since Weingartner & Draine (2001a) found that the photoelectric and 
photodesorption forces can be comparable to the direct radiation pressure force in some interstellar 
environments. 

Here, we estimate the photoelectric torque for two simple grain shapes, namely, prolate and 
oblate spheroids with axis ratios of 1.5. We also estimate the maximum magnitude for the pho- 
todesorption torque. We consider grains with silicate composition and a e g = 0.2/xm (a c g is the 
radius of a sphere with equal volume); such grains are known to be well-aligned in the diffuse ISM 
(e.g., Kim & Martin 1995). 

Photoelectrons typically originate within ~ 10 A of the grain surface and the ejection rate is 
proportional to the electric intensity |E| 2 . Thus, the field inside the grain must be determined to 
within high resolution (10 A out of 0.2/im) in order to compute the photoelectric torque, unless 
the near-surface field remains fairly constant on larger length scales. Draine & Weingartner (1996, 
1997) computed radiative torques using the discrete dipole approximation, in which the grain is 
approximated as a set of polarizable points. The size of the dipole array needed to achieve a 
resolution of 10 A out of 0.2^m is prohibitive. Thus, we adopt the point matching method (PMM), 
which can achieve high resolution at relatively low computational cost. In this approach, the 
internal and scattered fields are expanded in vector spherical wave functions and the expansion 
coefficients are found by imposing boundary conditions at a set of discrete points on the grain 
surface. 

The PMM simplifies considerably for grains with azimuthal symmetry; this motivated our 
choice of spheroids. However, these grains are too symmetric for radiation torques (including the 
recoil torques) to yield alignment, if the surface properties are uniform. The purpose of this paper 
is simply to estimate the magnitude of the recoil torques relative to the direct radiative torque, in 
order to determine whether or not recoil torques need to be included in future studies of radiation- 
driven alignment. Detailed scattering computations are required even for this simple goal, since the 
wavelength of the incident radiation is comparable to the grain size. In this case, the internal electric 
intensity cannot be simply estimated using geometrical optics and an attenuation coefficient. 

In £j2j we describe our implementation of the PMM, which is based on the treatments of 
Morrison & Cross (1974) and Al-Rizzo & Tranquilla (1995). $3] describes the computation of 
efficiency factors for extinction, scattering, radiation pressure, and radiative torque. Readers who 



- 3 - 



are interested only in the results should skip §£}2] and 13.11 In ££0 we describe our simple model 
for the photoelectric force and torque and present computational results. §5] briefly examines the 
photodesorption torque. We conclude ($6|) with a discussion of the implications of these estimates 
for future grain alignment studies. 



2. Point Matching Method 

In spherical coordinates, solutions of the vector Helmholtz equation are given by M nm (r, 9, 4>) 
and N nm (r,0,4>), with 
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where P(cos9) = P n n \cos9), P'(cos9) = dP^ (cos 0)/d(cos 9), k is the wave number {k = ko = 
2tt/X outside the grain and k = n re f&o inside the grain, where n re f is the grain index of refraction 
and A is the wavelength of the incident radiation), Z(kr) = j n (kr) (spherical Bessel function of 
the first kind) inside the grain and Z(kr) = hn (kyr) (spherical Hankel function of the first kind) 
outside the grain, and Z'{kr) = dZ (kr) / d(kr) . We adopt Jackson's (1999) convention for associated 
Legendre polynomials P™(cos#). 



'-nm 



The scattered electric (E sca ) and magnetic (H sca ) fields are given by 

co n 

E sca = ^ [ Q nm(l)NW +CW (2)M: 

n=l m=—n 

oo n 

H sca = -^ [onm(l)MW +a nm (2)N 

n=l m=—n 

and the internal fields are given by 
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where superscripts (s) and (w) refer to evaluation in the grain exterior and interior, respectively. 

We will take the z-axis as the symmetry axis of the spheroidal grain and assume that the 
incident radiation field wave vector k lies in the x-z plane with k = ko (cos 0qz — smO^x). We 
also adopt unit amplitude for the incident electric field E; nc . When Ei nc lies in the x-z plane (TM 
mode), 

Ei n c,TM = [f(cos#o srn 0cos<^+sin#ocos#)+#(cos#ocos#cos(/>— sin # sin #)—</> cos # sin (/>] exp(ik -r) , 

(11) 

H; nC) tm = [f sin sin (f> + cos sin <f> + <fi cos cf>] exp(iko • r) (12) 

with 

ko • r = /cor(cos #o cos — sin6*o sin 6* cos (ft) . (13) 
For the orthogonal polarization state (TE mode), 

Ei n c,TE = Hi nC) TM H inCj TE = ~ E; nCi TM • (14) 

At the surface, the tangential components of E and H and the normal components of D = eE 
(e is the dielectric function) and H are continuous. 

For the grain surface, the radius r as a function of polar angle is given by 

r{0) = a [cos 2 6> + C 2 sin 2 1/2 . (15) 

For prolate grains, a is the semimajor axis length and ( > 1 is the ratio of the semimajor axis 
length to the semiminor axis length. For oblate grains, a is the semiminor axis length and ( < 1 is 
the ratio of the semiminor axis length to semimajor axis length. Given the radius a e g of a sphere 
of equal volume, a = a e g £ 2//3 . 

In this case, the boundary conditions on the tangential field components become 

Einc,<p + -E-sca,^ = E- mt>( f, (16) 
H\nc,<f> + Hsc&,4> = Hint,<p (17) 

(1 - C 2 ) sin cos (E inCtr + E sca>r - E int , r ) + (cos 2 + ( 2 sin 2 0) (E- mcfi + E sca>0 - E intfi ) = (18) 
(1 - C 2 ) sin cos 6 (H iQCtr + H sca>r - H^ r ) + (cos 2 + Q 2 sin 2 9) (H inCt0 + H sc&fi - H int>e ) = . (19) 
The boundary conditions on the normal field components are 

(COS 2 + C 2 Sin 2 0) (Anc,r + Aca.r " Ant,r) + (( 2 " 1) sin COS 9 (Anc,0 + Aca,e " Ant,e) = (20) 

(cos 2 9 + C 2 sin 2 9) {H hlC)T + H sca<r - H int>r ) + (C 2 - 1) sin 9 cos 9 {H inc>e + H scafi - H hitfi ) = . (21) 
The incident radiation field can be expressed as follows: 

oo 

E r ,TM = exp(ifcor cos #o cos 9) i m [J m (u) sin 9q cos 9 — iJ' m (u) cos sin9] exp(im0) (22) 

m=— oo 
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E$,tm = — exp(z/cor cos 9q cos 9) i m [J m {u) sin Bq sin 6 + iJ' m (u) cos 6q cos 6] exp(im0) (23) 

m=— oo 

oo 

Eip^TM = exp(z/cor cos Oq cos 0) cos 6q i m mu~ l J m (u) exp(im0) (24) 

77l= — OO 

oo 

ff r ,TM = — sin6> exp(i/co^ cos #o cos 9) i m mvT x J m (u) exp(im<p) (25) 

m=— oo 

oo 

He,TM = — cos 6* exp(zfcor cos Bq cos 0) i m mu~ l J m {u) exp(im0) (26) 



m=— oo 
oo 



ff^TM = — exp(zfcor cos #o cos 0) i m+1 J' m {u) exp(im0) (27) 

m=— oo 

with u = — fcor sin^o sin6*, J m (u) the Bessel function, and J' m {u) = dJ m {u) / du. 

To find the expansion coefficients a nm , the field expressions in equations ([7]) through (|10p 
and ([22]) through ([27]) are substituted into the boundary conditions (fT6|) through (fl"9|) . In the 
resulting equations, each term contains a factor exp(im0) and no other 0-dependence. Multiplying 
by exp(— im'4>) and integrating over 0, we find the following equations, which much be satisfied for 
each value of m, 6, and i = 1 through 4: 

oo 4 

),C nm (i,j) ■ a nm (j) = D m (i) (28) 

n=|m| j=l 



with 
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C nm (l, 3) = -(1 - C 2 ) sin^ cos 0JV<$ r (0 = 0) - (cos 2 + C 2 sin 2 0)^/0 = 0) (31) 
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C nrn (3, 4) = n ie {C nm (l, 3) 
Cnm(4, 1) = C nm (2,2) 

Cnm{^i 2) = C nm (2, 1) 

Cnm(4, 3) = n re f C nrrt (2, 4) 
C„ m (4, 4) = n re fC nm (2, 3) 



(40) 
(41) 
(42) 
(43) 
(44) 



and 




(45) 
(46) 
(47) 
(48) 



In practice, we impose an upper cutoff for \m\, m max , and for n, n max . For each m, we find 
the expansion coefficients a nm by performing a least-squares minimization of 



We generally take 2n max to 3n max values of 9, evenly spaced in cos 9. 

The design matrix for the least-squares problem is a 4ng by 4(n max — \m\ + 1) matrix, made 
up from C nm (i,j; 9), and does not depend on 0q. We first apply a QR-decomposition to the design 
matrix. The computational time required to then complete the solution is small compared with the 
time to accomplish the factorization. Thus, several additional values of 9q can be treated without 
substantially increasing the computational time. The matrix equations often exhibit serious ill- 
conditioning. We work around this problem by implementing the solution using Mathematica, 
which supports arbitrary-precision arithmetic. 

Note that, when m changes sign, C(l,2), (7(1,4), C(2, 1), (7(2,3), (7(3,1), (7(3,3), C(4,2), 
(7(4,4), D (2), and -D(3) change sign, while the other components retain their sign. Thus, for the 
TM mode, 

a n -m(l) = a n , m (l) ; a„ _ m (2) = -a n , m (2) ; a n _ m (3) = a n , m (3) ; a„ _ m (4) = -a n>m (4) (TM) 



and for the TE mode, 

«n,-m(l) = — On,m(l) ! an,-m(2) = On,m(2) ! dn,-m(3) = — On,m(3) 5 a ni _ m (4) = a„ jm (4) (TE) 



As a result, it is not necessary to separately compute the expansion coefficients for m < 0. When 
m = 0, a„ jm (2) and a„ im (4) equal zero for the TM mode and a n>m (l) and a n ,m(3) equal zero for 
the TE mode. 




(49) 



(50) 



(51) 
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When 9q = 0, all of the azimuthal modes vanish, except for those with m = ±1. In this case, 
the above expressions for the incident field components (eqs. [22TI27P contain the ill-defined terms 
J' m (u) and J m (u)/u with u = 0. Equations (|45p through (|48p must be replaced with 

D^l) = D 1 (3) = --cos0exp(zfc o rcos0) (52) 

D x (2) = Di (4) = ~ exp(ik r cos 6») . (53) 

To verify the solution, we check that the boundary conditions (fT6|) through (f2Tj) are indeed 
satisfied. This is particularly powerful for conditions (|20|) and (f2Tj) . since these are not used in the 
solution. In addition, we check that the field intensity exhibits the following symmetry: |E| 2 (27r — 
(/)) = |E| 2 (0), and that the intensity is independent of 4> when 6>o = 0. 

Figures [1] through 0] show the internal and external (incident plus scattered) fields at the grain 
surface, for a couple illustrative cases. 

3. Efficiency Factors 

In this section, we describe the computation of efficiency factors for scattering, extinction, 
force, and torque. The calculation of the direct radiative torque is needed for comparison with the 
photoelectric and photodesorption torques. The other efficiency factors will be used to validate our 
PMM code ( O . 

3.1. Definitions and Formulae 

The scattering efficiency factor Q sca , is defined by 

-Psca = -FincQscaVraefj (54) 

where P sca is the scattered power and Fi nc is the incident flux. The scattered power is found by 
integrating the flux in the scattered fields over a sphere at infinity, yielding (Morrison & Cross 
1974; Al-Rizzo & Tranquilla 1995) 

. m>0 / t \ / \ i 

„ 4 v— v n[n + 1) (n + m) r . , , l9 , , , l9n . , 

2n +!(„-„)! [K-WI' + l—On i (») 

e m = 1 (2) for m = (m > 0). 

The extinction efficiency factor Q ext is defined in the same way as Q sca (eq. 154]) . except that 
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P sca is replaced by P a b s + -Psca, where P a b s is the absorbed power. From the optical theorem, 
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— Tm i 

(^OOcff 
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(56) 
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(57) 



(Al-Rizzo & Tranquilla 1995). The absorption efficiency factor Qabs = Qcxt — Qsca- 
The force, or pressure, efficiency factor Q pr is defined by 

Frad = TTOeff^radQpr (58) 

where F ra( j is the force on the grain and u ra( j is the incident radiation field energy density. The 
torque efficiency factor Qr is defined by 

Trad = TTOeff^radTT-Qr (59) 

where r ra d is the torque and A is the wavelength of the radiation. 

Farsund & Felderhof (1996) derived expressions for the force and torque when the incident 
and scattered waves are expanded in vector spherical wave functions. Since their conventions differ 
somewhat from ours, we briefly review some relevant results. First, the incident plane wave is given 
by 

Fine — ^ ^ \Pnm^nm ~\~ Qnm^^-nm] (60) 

n,m 

For the TM mode, 

Pnm = -d nm sin 9 P' (cos 6 ) ; q nm = d nm mP (cos 9 )/ sin 6 (61) 
and for the TE mode, 

Pnm = -imd nm P (cos 9 )/ sin 6 ; q nm = id nm sin 9 P' (cos 9 ) ; (62) 

^ = i"+'(-ir ftni"l1l ■ < 63 > 

n(n + 1) (n + my. 

For the TM mode, the coefficients that appear in eqs. (7.25) through (7.30) in Farsund and 
Felderhof are given by 

j _ -imd nm P(cos6 ) 
i -d nm sin6> .P'(cos6>o) 
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For the TE mode, 



"hum 
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(68) 



and the expressions for c^ nm and Cg nm are the same as for the TM mode, except for an additional 
negative sign when m < 0. 



To compute the efficiency factors Q pr and Qr, the expressions in equations (fBTj) through ([68]) 
are substituted into the Farsund & Felderhof (1996) expressions for the force and torque (their 
eqs. 7.25 through 7.30). 



3.2. Results 

In Figures [5] through [T0| we display computational results for the efficiency factors, for a 
prolate grain, with £ = 3/2, and an oblate grain, with £ = 2/3. We take a er T = 0.2/um and several 
values of 9q, the angle between the incident radiation field and the grain's symmetry axis. We adopt 
dielectric functions for "astronomical silicate" from Draine (2003). Here, and throughout the paper, 
we assume unpolarized incident radiation; thus, we average over the TE and TM modes. Although 
the force has a component perpendicular to ko, we only display Q pr = Q pr • ko. The torque always 
lies along y, and vanishes when cos#o = or 1. (Recall the definition of the coordinate system, in 
the paragraph following eq. [TUl ) 

For various values of £, A, and we have also computed the efficiency factors using the discrete 
dipole approximation code DDSCAT (version 6.1, Draine & Flatau 2004). In all cases, the results 
were identical, to within uncertainties associated with incomplete convergence. We encountered 
no obstacles in convergence with the PMM code, though in some cases we had to take n max = 80 
to converge to within 1%. For some wavelengths, we were unable to converge to within 10% with 
DDSCAT, due to limitations on the size of the dipole arrays when using the GNU Fortran compiler. 
This successful comparison, using two completely different methods, validates both our PMM code 
and DDSCAT. 

For a given incident radiation field, the spectrally averaged force and torque efficiency factors 
are given by 

Qpr = / Qpr u v dv (69) 

«rad J 

Qr = y~ — / QrXu v du , (70) 

A «rad J 
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where v is the frequency, u v is the specific energy density, and A = / \ u v dv /tt ra d- In Figure fTTj we 
display Qr,y for the average interstellar radiation field (ISRF) spectrum in the solar neighborhood, 
as estimated by Mezger, Mathis, & Panagia (1982) and Mathis, Mezger, & Panagia (1983). (See 
eq. 31 in Weingartner & Draine 2001b for a convenient representation of the ISRF.) Interestingly, 
the torque due to radiation with A > 1/mi substantially cancels the torque due to radiation with 
A < 1/im. For the ISRF, u rad = 8.63 x 10~ 13 erg cm" 3 and A = 1.20/xm. 

Both the wavelength-dependent and spectrally-averaged efficiency factors are comparable in 
magnitude to those for the irregular grain studied by Draine & Weingartner (1996); compare figs. [9] 
and [10] with their fig. 7 and fig. [TT] with their fig. 12. 

For a spherical silicate grain with a e s = 0.2/im exposed to the ISRF, Q pr = 0.82. For spheroids 
with £ = 2/3 (3/2), Q pr ranges from 0.78 to 0.89 (0.73 to 0.85), with larger values corresponding 
to larger cross-sectional areas. The transverse component of Q pr takes values as high as ~ 0.1. 



4. Photoelectric Forces and Torques 

4.1. A Simple Model 

For simplicity, we will assume that (1) the grains are electrically neutral, (2) the photoelectrons 
emerge along the surface normal n, and (3) the photoelectron energy E e = hv — W, where hu is 
the incident photon energy and W is the work function of the grain material. Weingartner & 
Draine (2001b) estimated that silicate grains with a e ff = 0.2//m in the cold neutral medium have 
a slight positive charge. The rate at which photoelectrons are ejected and the average energy per 
photoelectron are both lower for positively charged grains than for neutral grains. Thus, each of 
these simplifications results in an overestimate of the force and torqueQ This is appropriate for 
our exploratory study, since we are only attempting to ascertain whether or not the photoelectric 
torque may be important enough, compared with the direct radiative torque, to warrant further 
consideration in studies of grain alignment. We will adopt W = 8eV and the photoelectric yield 
Y (i.e., the probability that a photoelectron is ejected following the absorption of a photon) 

these are the estimates of Weingartner & Draine (2001b) for bulk silicate. On theoretical grounds, 
it is expected that the yield and ionization potential for a grain with a e g = 0.2/im deviate little from 
their bulk values (e.g., Watson 1972, 1973). However, scant experimental evidence on photoelectric 
emission from sub-micron grains is available. The recent experiment of Abbas et al. (2006) found 
that the yields of such grains substantially exceed the bulk yield. The work function and yield have 



1 In the rare environments where the grains are highly charged, these may be overestimated by over an order of 
magnitude (Weingartner & Draine 2001a). 
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not been well characterized experimentally, even for bulk silicate. Thus, our results for the force 
and torque will only be rough estimates, and we are justified in making the above simplifications. 

Following Kerker & Wang (1982), we will take the photoemission rate as a function of position 
on the surface proportional to the internal field intensity evaluated at the grain surface, |E™ t rf | 2 . 
With this assumption, as well as those in the preceding paragraph, the photoelectric force F pe and 
torque r pe are given by 

fpe — TTQ-eft ^radQpr, pe 



~pe — TTflefr^rad-— Qr, pe 
Z7T 



with 



Qpr, pe 
Qr,pe = 
A — 



cQabs 
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I l E int 



surf|2 nd5 



J\B^\ 2 dS 

J\EZ'\ 2 rxhdS 
r a eS J\E^dS ' 

p e is the photoelectron momentum. 
The unit surface normal 

n = [1 + (C 4 - l)sin 2 6*]~ 1/2 [C 2 sin 9 (x cos <p + y sin 4>) + zcos9] 

and the area element 

lV2 



dS 



r 2 sin l 



1 + 



;i-c 2 )- 



sin 



cos 



d0< 



(72) 
(73) 

(74) 
(75) 
(76) 

(77) 



(78) 



(79) 



4.2. Results 

The photoelectric torque, like the direct radiative torque, always lies along y and vanishes 
when cos#o = and 1. Figure [12] displays the y-component of the torque asymmetry factor Ar for 
two values of the incident radiation wavelength A, both lying within the relevant 8eV to 13.6 eV 
spectral range. 

In estimating the maximum value of the photoelectric torque, we will take |^4r,y| ~ 0.05, 
independent of A (fig. [T2|) . We will also take Q a bs ~ 1, a reasonable approximation when the 
photon energy exceeds 8eV (figs. [5] and [6]). Defining the spectrally-averaged efficiency factor as 
in equation (|70|) . we find that 0r,pe, y ~ 3.8. The energy density in the ISRF between 8eV and 
13.6 eV is u(> 8eV) = 3.86 x 10 _14 ergs _1 and A = 0.124/im for this spectral range. From Figure 



-12- 



[TT| |Qr,y| reaches values as high as ~ 0.05 for the ISRF. Comparing the radiative and photoelectric 
torques at the Q$ for which they peak, we find that r po ~ 0.35 r ra d. 

If the photoelectric yield varies across the grain surface, then the torque is no longer confined 
to lie along y, raising the possibility that even symmetric shapes like spheroids could experience 
radiation-driven alignment (Purcell 1979). In this case, equation (I77p for Ap is modified, with the 
yield Y appearing in both integrals. The z-component of Ap remains zero, since the z-component 
ofrxn always vanishes. Specifically, we situated a single large spot on the grain surface where 
the yield is enhanced by 10%. When the spot is on the illuminated side of the grain, Ar, y changes 
by as much as ~ 20% and \Ap x \ (and |-Ar )3/ | for cos#o = 0) 1) reach values as high as ~ 0.005. 
Thus, for real interstellar grains, the aligning torque associated with non-symmetric shape likely 
dominates any aligning torque associated with non-uniform photoelectric yield. 

The forward-direction force asymmetry parameter A pT = A pr • ~ 0.4 to 0.6, while that for 
the transverse direction ranges from zero to ~ 0.15. With A pr ~ 0.5 and Qabs ra 1> the ISRF 
spectrally-averaged efficiency factor Q pr , P e ~ 3.7, yielding F pe 0.2F rax ±. 

5. Photodesorption 

Although uncertainties abound in modelling the photoelectric force and torque, photodesorp- 
tion presents even greater challenges, as described in §4 of Weingartner & Draine (2001a). A major 
question is whether or not adsorbed H atoms can diffuse across the surface, by either thermal bar- 
rier hopping or quantum mechanical tunneling. This depends on the poorly known surface-adatom 
binding energy. 

For uniform surface properties, the magnitude of the photodesorption torque increases with 
increasing surface coverage of adatoms. Since we seek to estimate the maximum plausible photodes- 
orption torque, we assume complete coverage across the entire surface. This results, for example, 
if the adatoms do not desorb and the rate at which gas-phase atoms collide with the grain greatly 
exceeds the removal rates (due to H2 formation as well as photodesorption). 

Photodesorbed H atoms are produced at a rate < R® d S/l 2 , where R® d is the photodesorption 
rate per adatom, S is the grain surface area, and I 2 is the surface area per binding site. We adopt 
the crude estimate from Weingartner & Draine (2001a) of R® d ~ 2 x 10~ 10 s _1 for grains exposed to 

the ISRF, and we take I 2 ~ 10 A 2 . We also assume that the kinetic energy of photodesorbed atoms 
is comparable to that of photoelectrons, i.e., ~ eV . Finally, we assume that the photodesorption 
rate as a function of position on the surface is proportional to the electric intensity just above the 
surface. 

The ratio of the photodesorption torque F pc j to the photoelectric torque can then be simply 
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estimated as 




S R%hu 



( 




1/2 



|Ar(pd)| 
|A r (pe)| 



(80) 



where S « 47ra 2 ff for (" = 3/2 and 2/3, mjj is the proton mass, Ap(pe) is the asymmetry parameter 
given in equation (|77p . and Ar(pd) is identical, except that the internal electric field E™ t is 
replaced with the external electric field at the surface, E|" f . For A = 0.15//m, -Ar )3/ (pd) differs very 
little from Ar,y(pe), but for A = 0.1/xm, the external asymmetry factor is as much as three times 
smaller than the internal asymmetry factor. With Qabs ~ 1> Y ~ 0.07, and ^r(pd) ~ Ar(pe), 
Tpd 0.3 r pe . 

The external and internal pressure asymmetry factors are also generally comparable, though 
the former can be as much as a factor of 2 smaller than the latter. Thus, for the models considered 
here, F p( ^ and F pe are of comparable magnitude, as Weingartner & Draine (2001a) found for spheres. 



The goal of this work was to determine whether or not recoil torques must be included in 
investigations of radiation-driven grain alignment. Our model for the photoelectric torque neglected 
a few effects that would tend to suppress the torque ( §4.ip . We found r po « 0.35r ra( j, suggesting 
that the photoelectric torque may actually contribute at about the 10% level compared with the 
direct radiative torque. 

However, the photoelectric yield for sub-micron silicate grains is highly uncertain. If the recent 
experimental results of Abbas et al. (2006) are accurate, then the yield may be substantially larger 
than we have assumed. This would increase the torque, but not in direct proportion, since higher 
yield implies higher positive grain charge, which partially suppresses the torque. Yet, we cannot 
rule out the possibility that the photoelectric torque is of the same order of magnitude as the 
radiative torque for spheroids. Clearly, additional experiments are needed to better characterize 
photoelectric emission from sub-micron silicate grains. 

We must also acknowledge the possibility that less symmetric grains may experience larger 
recoil torques than spheroids, even though the direct radiative torques for the spheroids examined 
here are comparable in magnitude to that for the irregular grain studied by Draine & Weingartner 
(1996). The computational demands associated with the PMM are much more severe for grains 
lacking azimuthal symmetry. However, we have found that the near-surface fields for spheroids are 
extremely smooth (figs. CD through 0]), suggesting that the discrete dipole approximation may yield 
sufficient resolution for accurate evaluation of the recoil torques. We will pursue this possibility in 
future work. 

The photodesorption torque appears to contribute at most at the 10% level, although it may be 
more important if our adopted value for the photodesorption rate, R pd , is too small. In making our 
simple estimate, we assumed that the surface-adatom binding energy is relatively large (> leV), 



6. Conclusion 
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but smaller binding energies (as suggested by, e.g., the recent experiment of Perets et al. 2007) 
may result in a smaller torque. Futher experimental work characterizing the surfaces of amorphous 
silicates, and the photodesorption rates, is needed before the torque can be reliably estimated. Thus, 
for now, we recommend that both recoil torques be omitted from detailed studies of radiation-driven 
alignment. 

We are grateful to Bruce Draine and an anonymous referee for helpful comments. JCW is 
a Cottrell Scholar of Research Corporation. Support for this work, part of the Spitzer Space 
Telescope Theoretical Research Program, was provided by NASA through a contract issued by the 
Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. 
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Fig. 1. — The intensity I = |E| 2 of the internal field (averaged over the TE and TM modes and 
normalized to the incident intensity Jo) at the grain surface, for an oblate (£ = 2/3) silicate grain 
with a c ff = 0.2/xm, incident wavelength A = 0.1/xm, and angle between the incident radiation and 
the grain symmetry axis 9q = cos" 1 0.6. 
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Fig. 2. — Same as fig. HJ except that the external field intensity at the grain surface is plotted. 
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Fig. 3. — Same as fig. HJ except that A = 1/xm. 
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Fig. 4. — Same as fig. EJ except that A = 1/xm. 
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Fig. 5. — Qabs for prolate (( = 3/2) grains with a e g = 0.2/xm and various values of 9q. 
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Fig. 6. — Q a bs for oblate (( = 2/3) grains with a e ff = 0.2/um. 
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Fig. 7. — Q pr for prolate (( = 3/2) grains with a e g = 0.2/j.m. 
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Fig. 8. — Q pr for oblate (£ = 2/3) grains with a e fr = 0.2//m. 
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Fig. 9. — XQr,y for prolate (( = 3/2) grains with a c g = 0.2//m. 
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Fig. 10. — XQr,y for oblate (( = 2/3) grains with a c g = 0.2//m. 
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Fig. 11. — The torque efficiency factor averaged over the MMP radiation field, Qr,y, for prolate 
(C = 3/2) and oblate (C = 2/3) grains with a e s = 0.2/xm. 
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Fig. 12 — 



The y-component of the torque asymmetry factor Ap. 



